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BLIND ADAPTIVE FILTERING METHOD FOR RECEIVERS OF 
COMMUNICATION SYSTEMS 

FIELD OF THE INVENTION 
The present invention relates to a method of blind adaptive filtering for receivers 
of communication systems. 

BACKGROUND OF INVENTION 

The present invention relates to communications equipment, and more 
particularly to a blind adaptive filtering algorithm for receivers of communication 
systems without need of training sequences. 

In wireless multi-user communication systems, antenna array based non-blind 
linear minimum mean square error (LMMSE) receivers utilize long training sequences 
for the suppression of multiple access interference (MAI) due to multiple users sharing 
the same channel and the suppression of intersymbol interference (I SI) due to channel 
distortion in wired communications and multi-path propagation in wireless 
communications. The MAI suppression and ISI suppression are crucial to the 
performance improvement of wireless communication systems. However the use of 
training sequences not only costs system resources but also involves synchronization 
issues and other cumbersome procedures. Therefore, receivers for blind mode (without 
need of training sequences) with good performance (of suppression of both MAI and 
ISI), and low complexity are of extreme importance in communication systems. 

In digital communication systems, by sampling or matched filtering, the received 
continuous-time signal can be transformed to a discrete-time signal which can be 
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approximated as the output signal of a discrete-time multi-input multi-output (MIMO) 
linear time-invariant (LTI) system on which a variety of detection and estimation 
algorithms are based. The proposed blind adaptive filtering algorithm for receivers of 
communication systems is also based on the discrete-time MIMO LTI system model. 

Blind equalization (deconvolution) of an MIMO LTI system, denoted H[n] (P x 
K matrix) is a problem of estimating the vector input u[n] = (u } [n], u^n], u K [n]) T 
(K inputs), where the superscript 'T' denotes the transposition, with only a set of non- 
Gaussian vector output measurements x[n] = (^[n], x^n], Xp[n]) T (P outputs) as 
follows and as shown in Figure 1 : 



where '*' denotes the discrete-time convolution operator, hy[n] is the jth column of 
H[n], and w[n]=(if 1 [n.],if 2 [n.],. . .,Wp[n]) T (P x 1 vector) is additive noise. For ease of 
later use, let us define the following notations 

cum{y l ,y 2 , ,y P } pth-order cumulant of random variables y x ,y 2 ,- -,yp 

cum{y :/?,...} = cum^ = y, y 2 = y, . ., y P = y, . ■■ } 

C pq {y) = cum{y : p, y: q} {y is complex conjugate of y) 

v j =(v j [L l lv J [L l +ll. .,Vj[L 2 ]f ((L = 1^-^+1) x 1 vector) 



x[n] = H[n] * u[n] + w[n] 

= J H[fc]u[n - fc] + w[n] = f h^n] * u^] + w[n] 



0) 




v = 



Xj [n] = {Xj [n - Lj ], Xj [n - L, - 1], . . ., x j [n - L 2 ]) T 



x[n] = (x 1 r [n],x 2 r [n] ) ... ) x;[n]) r 



R =E[x[n]x T [n]] 



(expected value of x* [n]x T [n]) 
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Assume that we are given a set of measurements x[n], n = 0, 1, . . ., N-l modeled 
by (1) with the following assumptions: 

(Al) Uj[n] is zero-mean, independent identically distributed (i.i d.) non-Gaussian with 
variance E[\ Uj [n] | 2 ] = and (p+ q)th-order cumulant C p q {uj [n] } ^ 0 , 
and statistically independent of u k [n] for all k ^ j , where p and q are 
nonnegative integers and p + q > 3 . 

(A2) The MIMO system H[n] is exponentially stable. 

(A3) The noise w[n] is zero-mean Gaussian and statistically independent of u[n]. 

Let v[n] = (v } [n], v 2 [ri\, . ., v P [n]) T denote aP x l linear finite impulse response 
(FIR) equalizer of length L = L 2 — L } + 1 for which v[n] * 0 (P x 1 zero vector) for n = 
L b Lj+1, . . L 2 . Then the output of e[n] of the FIR equalizer v[n] can be expressed as 

e[n] = j^vj[n] * x^n] = f>J*,[n] =v r x[n] (2) 

Tugnait (see "Identification and deconvolution of multichannel linear nonGaussian 
processes using higher-order statistics and inverse filter criteria," IEEE Trans. Signal 
Processing, vol. 45, no. 3, pp. 658-672, March 1997) proposed inverse filter criteria 
(IFC) for blind deconvolution of MIMO systems using second- and third-order 
cumulants or second- and fourth-order cumulants of inverse filter (i.e., equalizer) output 
Chi and Chen (see "MIMO inverse filter criteria and blind maximum ratio combining 
using HOS for equalization of DS/CDMA systems in multi-path," Proc. Third IEEE 
Workshop on Signal Processing Advances in Wireless Communications, Taoyuan, 
Taiwan, March 20-23, 2001, pp. 114-1 17, see also "Cumulant-based inverse filter criteria 
for MIMO blind deconvolution: properties, algorithms, and application to DS/CDMA 
systems in multi-path," IEEE Trans. Signal Processing, vol. 49, no. 7, pp. 1282-1299, 
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July 2001, and "Blind beamforming and maximum ratio combining by kutosis 
maximization for source separation in multi-path," Proc. Thrid IEEE Workshop on 
Signal Processing Advances in Wireless Communications, Taoyuan, Taiwan, March 20- 
23, 2001, pp. 243-246) find the optimum v by maximizing the following IFC 



where p and q are nonnegative integers and p + q > 3 , through using gradient type 
iterative optimization algorithms because all J p>q (v ) are highly nonlinear functions of 
v (without closed-form solutions for the optimum v ). Note that the IFC given by (3) 
include Tugnait's IFC for (p,q) = (2,1) and (p,q) = (2,2) as special cases. Under the 
following conditions: 

(al) signal-to-noise ratio (SNR) is infinity, 

(a2) the length of equalizer is infinite, 
the optimum equalizer output turns out to be one of the K input signals except for an 
unknown scale factor and an unknown time delay (i.e., the optimum equalizer is a 
perfect equalizer for one of the K input signals). 

Yeung and Yau (see "A cumulant-based super-exponential algorithm for blind 
deconvolution of multi-input multi-output systems" Signal Processing, vol. 67, no. 2, pp 
141-162, 1998) and Inouye and Tanebe (see "Super-exponential algorithms for 
mulichannel blind deconvolution," IEEE Trans. Signal Processing, vol. 48, no.3, pp. 
881-888, March 2000) also proposed an iterative super-exponential algorithm (SEA) for 
blind deconvolution of MIMO systems. The iterative SEA updates v at the Ah iteration 

by 
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y,= °i (4) 

where || a || denotes the Euclidean norm of vector a and 

=cum{e H [n]:r,(e H [n])* :s-l,x*[n]} (5) 
in which r and s— 1 are nonnegative integers, r + s>3 and e^j/z] is the equalizer 
output obtained at the (i — l)th iteration. Again, under the conditions (al) and (a2), the 
designed equalizer by the SEA is also a perfect equalizer (for one of the K input signals) 
with a super-exponential convergence rate. A remark regarding the distinctions of 
gradient type IFC algorithms and the SEA are as follows: 

(Rl) Although the computationally efficient SEA converges fast (with a super- 
exponential convergence rate) for SNR =oo and sufficiently large N, it may 
diverge for finite N and finite SNR. Moreover, with larger computational load 
than updating the equalizer coefficients v by (4) at each iteration, gradient 
type iterative IFC algorithms (such as Fletcher-Powell algorithm (see D.M. 
Burley, Studies in Optimization, Falsted Press, New York, 1974)) always 
converge slower than the iterative SEA for p+q= r+s as x[n] is real and for 
(p,q) = (r,s) as x[n] is complex. 
Therefore a need exists for a fast IFC based algorithm with performance similar 

to that of SEA with guaranteed convergence for finite SNR and finite data length, and is 

therefore suitable for a variety of practical uses. 



SUMMARY OF THE INVENTION 
A method according to the present invention is proposed for fast iterative MIMO 
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blind deconvolution in a communications receiver using cumulant based inverse filter 
criteria with super-exponential convergence rate and low computational load. At each 
iteration the method updates the equalizer coefficients v using the same update 
equation as used by the super-exponential algorithm. If the value of the associated 
inverse filter criteria (IFC) is larger than that obtained at the immediately previous 
iteration, then the method goes to the next iteration. Otherwise the equalizer coefficients 
v are updated using gradient type optimization algorithms such that the value of IFC 
increases. As the algorithm converges, the optimum equalizer v[n] is obtained and its 
output e[n] is an estimate of one of the K driving inputs of the MIMO system. The 
convergence can be guaranteed and the number of iterations for convergence is usually 
between two and ten. 

These and other features, which characterize the invention, are set forth in the 
claims annexed hereto and forming a further part hereof. However, for a better 
understanding of the invention, and of the advantages and objectives attained through its 
use, reference should be made to the drawings, and to the accompanying descriptive 
matter, in which there is described exemplary embodiments of the invention 

BRIEF DESCRIPTION OF DRAWINGS 
FIGURE 1 shows the multiple-input multiple-output (MIMO) system 
FIGURE 2 shows the method of signal processing at each iteration of Algorithm 1 
FIGURE 3 shows the method of signal processing of each stage of the MSC procedure 

DETAILED DESCRIPTION 

In the following detailed description of the preferred embodiments, reference is 
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made to the accompanying drawings which form a part hereof, and in which is shown 
by way of illustration specific preferred embodiments in which the invention may be 
practiced. The preferred embodiments are described in sufficient detail to enable these 
skilled in the art to practice the invention, and it is to be understood that other 
embodiments may be utilized and that logical, changes may be made without departing 
from the spirit and scope of the present invention. The following detailed description is, 
therefore, not to be taken in a limiting sense, and the scope of the present invention is 
defined only be the appended claims. 

A key property of the optimum equalizer v[n] using IFC, that can be proven true 
for any signal to noise ratio (SNR), is that the optimum v[n] with finite L using IFC and 
the one obtained by the SEA are the same (up to a scale factor) for p+ q = r+ s > 3 as 
x[n] is real and for p = q = r = s > 2 as x[n] is complex. Based on this property and 
(Rl), the fast iterative MIMO blind deconvolution algorithm using IFC, called 
Algorithm 1, which processes the given x[n], n=0, 1, . . ., iV— 1, to estimate one of the 
K driving inputs u,[n], u^n], . . , u^n] comprises two steps as follows: 

In the first step given v 7-1 and e^fn] obtained at the (I— l)th iteration, v 1 at 
the Ih iteration is obtained by the following two steps and as shown in Figure 2: 
v j is obtained by using the following SEA equation 



with r+s = p+q as x[n] is real and with r = s = p = q as x[n] is complex and obtain 
the associated eXn], n = 0, 1, N—l, where a>0 is a parameter such that 



■T-1 



(6) 



v 




E[\ ei [nf]=a 2 
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In the second step, if J p g (v 7 ) > J p q (y go to the next iteration, otherwise 
update v 7 through a gradient type optimization algorithm such that 
J P ,q ( v / ) > J p,q (y i-\ ) and obtain the associated e£n], n = 0, 1 , . . . , N— 1 . 

As x[n] is complex for p = q = r = s and x[n] is real for p-q = r+s, it can be 
shown that 



where d j_ { has been obtained in step one (see (5)) and R is the same at each iteration, 
indicating simple and straight-forward computation for obtaining the gradient 
dJ p q (v )/ 5v in step two. 

The method of Algorithm 1 is an iterative batch processing IFC algorithm with 
convergence speed, computational load similar to those of SEA due to step one of 
Algorithm 1. Moreover, its convergence can be guaranteed because J p , q (y j) (which 
is bounded) increases at each iteration due to step two of Algorithm 1. An initial 
condition v 0 (eg., the one associated with v i [n] = S[i- j]S[n-l], i = 1, 2, . , P 
where 1 < j < P,L t <l < L 2 and S[n] = 0 for n ^ 0 and £[0] - 1 ) is needed to 
initialize Algorithm 1 and the resultant equalizer output e[n] (after convergence) is the 
obtained estimate of one of the K driving inputs u^[n], u^n], u K [n] of the MEMO 
LTI system. The parameter a in equation (6) plays the role of restricting the dynamic 
range of the equalizer output to avoid overflow problem in DSP implementation. 

Algorithm 1 can be applied wherever the unknown PxK system H[n] 
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approximates an LTI system or a slowly linear time-varying system, and that Algorithm 
1 can be applied as H[n] approximates either a single-input multi-output (SIMO) LTI 
system (i.e., P>1,K= 1) or a single-input single-output (SISO) LTI system (i e., P = K 
= 1). Hence all the estimates u,[n],u 2 [2], ...,u K [n] can be obtained by Algorithm 1 
(possibly in a non-sequential order) through the multistage successive cancellation 
(MSC) procedure that includes the following two steps at each stage as shown in Figure 
3. First, find an input estimate, said u^n] (where I is unknown), using Algorithm 1, 
and then obtain the associated channel estimate by 

• E[x[n + fcKM 

£[|",M! 2 ] 

In the second step update x[n] by x[n ] - h l [n ] * u l [n ], i.e., cancel the 
contribution of u t [n] in x[n] 

In conjunction with the MSC procedure, Algorithm 1 can be applied wherever the 
given non-Gaussian measurements x[n] can be approximated as the outputs of an 
MLMO system (channel) given by (1) without the need of channel information and 
training sequences. The following are some possible applications of the above algorithm 
although the invention is by no means limited to these applications. 

A smart antenna processing means is used in wireless communications when a 
communications receiver has an antenna array. The smart antenna processing means 
provides an estimate of a signal received from a particular user. For the scenario of K 
user signals arriving at an antenna array of P elements, the received P x 1 signal vector 
x(i) can be modeled as 

K M k 

X ( f ) = E 2>/- a (^m K - ) + W(t) (9) 
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where 6 km , a km and are the angle-of-arrival (AO A), propagation loss and 
propagation delay associated with the mth path of user k, respectively; M k denotes the 
total number of paths of user k; a(9 km ) denotes the array steering vector associated 
with 9 km ; w(t) is the noise vector; and 

s*(0= Z u k [n]p k (t-nT) (10) 

is the baseband signal of user k, where u k [n], p k (f), and T are the symbol sequence 
(non-Gaussian such as binary sequence of +1 and — 1), signature waveform and 
symbol duration, respectively. In time division multiple access (TDMA) systems, 
usually p k (t) = p(f) (pulse shaping function) for all k. In direct sequence code division 
multiple access (DS/CDMA) communications systems, the signature waveform can be 
chosen as 

p k (t) = -^ZOc in¥(t -nT c ) (11) 

where P is the spreading gain, T= TIP is chip period, c k [n] is a binary pseudo random 
sequence of {+1, — 1 }, and <j>(t) is the rectangular chip pulse of magnitude equal to 
unity within the interval t e (0, T c ) . 

The received signal vector x(f ) given by equation (9) can also be expressed as 

x(t) = Z Zu fc [n]h fc (t-nT) + w(t) (12) 

where 

h(0= lla km a{e km ) Pk {t-T krn ) (13) 

can be called the symbol-level space-time channel impulse response for user k. In 
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most wireless communication environments, the channels h k (t), k = I, . ., K can be 
assumed to be linear and time-invariant during the period of interest 

In the case of TDMA, the estimation of u^n], u^n], u k [n] is based on the 
symbol waveform. Denote A T = T/Q T as the sampling period, where Q T >1 is the 
factor of over-sampling. Sampling x(f) at t = nT -iA T , n e (-oc,+oo) 3 equation (12) 
becomes 

K oo 

x r (nT-iA r ) = ^ ^«Jfi-d]h)c(dT -iA T ) + w T (nT -iA T ), i = 0, 1, ...,Q T - 1 

(14) 

and an MIMO system model as given by equation (1) can be easily obtained as 



K oo 



(15) 



= ZH[d]u[n-d] + w[n] 



(16) 



where 



x[n] 



x T (nT) 
x T (nT- A T ) 



(17) 



x r (nT-(Q r -l)A r 



J(pg r )xi 



\L k [d] = 



h k (dT) 
h k (dT-A T ) 



(18) 



h fc (dT-(Q r -l)A r )_ 



'J(pg r )xi 



11 



FILE: 7934USF.RTF 



w[n] = 



w r (nT) 
w T (nT - A T ) 



(19) 



w r (nT-(Q T -l)A T )_ 



'J(P(? r )xl 



and 



H[d] = [h 1 [d],h 2 [d],...,li K [d]] 



(20) 



is a (PQt) x -K" system. Therefore, all the user's symbol sequences u[n] = (u x [n], 
u^n], . .., u^_n]) T can be detected using Algorithm 1 through the MSC procedure 
followed by threshold decision. Therefore Algorithm 1 performs as a blind space-time 
multi-user detector as P > 2, whereas as P = /, Algorithm 1 performs as a blind multi- 
user detection algorithm (without space diversity). 

In the CDMA case, the estimation of u^[n], x^[n], .. , u^n] is based on the 
signature wave form and related matched filtering outputs. Denote A c = T c I Q c as the 
sampling period, where Q c > 1 is the factor of over- sampling. Sampling x(t) at 
t = lT c - iA c , i = 0, 1, Q c - 1, equation (12) becomes 



x c (lT c -iA c ) = £ f^u k [n]h k (lT c -iA c -nT) + w c (lT c -iA c ) (21) 



By taking x c (/T c -iA c ), i = 0, 1, Q c — 1, one can obtain the (PQ C ) x 1 signal 



vector 



ic ( /T c ) = Z 2X Milk ( ?T c - nT) + w c (IT C ) 



(22) 



where 



x c {IT C ) = [x T c (IT C ), x T c (lT c -A c \...,x T c (IT C - (Q c - l)A c )] T 



(23) 
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!i (IT C ) - [h£ (IT C Xh T k (lT c -A c ),...,h T k (IT C - (Q c - \)A c )f (24) 
w c (IT C ) = [ w T c (IT C ), w T c (IT C - A c ), . . w I (IT C - (Q c - 1) A c )] r (25) 

Assume that the signature waveform is chosen as equation (11). Then the discrete 

signature waveform of the /cth user is given by 

Pk (lT c ) = -j=Y.c k [nmiT c -nT c ) = ^=c k [l] (26) 

(since <f>(£) is the rectangular pulse) 
where 0 < I < P- 1 , and the output signal vector after despreading with the signature 
waveform of user k$ can be expressed as 

^'VT c ) = X*c(lT c -iT c )p ku (PT C -iT c ) (27) 
= E Z" fc [i]qJ° (Zr e -nT) + w? (ZT c ) (28) 



where 



q ;:° (iT c )=^fh k (iT c -iT c )c^\i] 

(ZT C ) = § w c (ZT C - zT c )c fcQ [i] (30) 

cJi] = cJP-/L 0<0<P-1 (31) 
Let IT C = mT - ijT c , j = 1, 2, . . ., J < P, be the J time indexes associated with the J 

largest values of Xm=o|^c ( mT ~ iT c )f > * = °> ! > • • •> P ~ 1 • Then a symbol-level 
MIMO signal model as given by (1) can be obtained as 

x /co [m] = 2 J u fc [n]q^(mr-nT) + w fc °[m] (32) 



(29) 
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]TQ[rn-n]u[n] + w fc ° [ m ] 



(33) 



and 



x ,c ° [m] = 



q;;°[mT] = 



°[m] = 



"^(mr-t-Tj" 

q £° K-i,r c )' 

q*° (mT-i 2 T c ) 
q*° (mT-i.TJ 

wZ?(mT-i 2 T c ) 

w k c °(mT-ijT c )^ 



(34) 



(35) 



(36) 



(37) 



Q[m] = [q; co (mT), q* 0 (mT), qj° (mT)] 
is a (JPQ C ) x X system. Therefore, all the user's symbol sequences u[n] = (u x [ri\, 
u^n], .. , u K [n]) T can be detected using Algorithm 1 through the MSC procedure 
followed by a threshold decision. 

Another application for the algorithm of the invention is for blind equalization in 
a communications system. The baseband received signal x(t) in the presence of M paths 
in wireless communications (or channel distortion in wired communications) can be 
expressed as 



x(t)=£ a m s(t-r m ) + w(t) 



(38) 
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where a m and r m are the amplitude and propagation delay associated with the mth 
path, respectively; w(t) is the noise; and 

s(t)= f] u[n]p(t-nT) (39) 

is the user's signal where u[n\, p(f) and Tare the symbol sequence (non-Gaussian such 
as binary sequence of +1 and — 1), signature waveform and symbol duration, 
respectively. 

By sampling x(t) at t = nT-i(T/Q), i = 0, 1, Q-l, where Q> 1 is the 
integer factor of over-sampling and then stacking the Q over-sampling vector samples, 
one can obtain an MIMO model given by (1) where H[n] is a Q x 1 system, and 
therefore, the user's symbol sequence u[n] can be obtained using Algorithm 1 followed 
by threshold decision. Therefore, Algorithm 1 performs as a blind equalizer for the 
suppression of either multi-path effects in wireless communications or channel 
distortion in wired communications. 

A further application of the invention is in blind RAKE receivers In the signal 
processing procedure of RAKE receivers, assume that there are P transmission paths 
(fingers) associated with the desired user and the P x 1 vector x[n] = (x,[n], x,[ft], 
x P [n]) T includes P estimates of the desired user's symbol sequence s[n] associated with 
the P respective transmission paths. The x[n] can be modeled as 

x[n] = as[n] + w[n] (40) 
where a is a P x 1 (unknown) column vector of P signal (complex) amplitudes of s[n] 
associated with the P fingers, and the P * 1 vector w[n] is noise. We desire to estimate 
s[n] from x[n] with no need of any information about a. 

Note that the P x 1 system given by (40) is a special case of the MIMO model 
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given by (1) with H[n]* 0 (P x 1 zero vector) only for n= 0. Algorithm 1, with L, = L 2 
= 0 (thus L = I, and v = v[0]), can be employed to obtain an estimate 

s[n] =v r x[n] = (v r a)s[n] +v r w[n] (41) 
that can be shown to have the maximum of 

SMt= a£f)*£i (42) 

E[\\ T w[n]\ 2 ] 

attained by the non-blind LMMSE estimator of s[n] that requires the vector a given in 
advance. 

Various additional modifications may be made to the illustrated embodiments 
without departing from the spirit and scope of the invention. Therefore, the invention 
lies in the claims hereinafter appended. 
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